Constraints on stable equilibria with fluctuation-induced forces 
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We examine whether fluctuation-induced forces can lead to stable levitation. First, we analyze a 
collection of classical objects at finite temperature that contain fixed and mobile charges, and show 
that any arrangement in space is unstable to small perturbations in position. This extends Earn- 
shaw's theorem for electrostatics by including thermal fluctuations of internal charges. Quantum 
fluctuations of the electromagnetic field are responsible for Casimir/van der Waals interactions. Ne- 
glecting permeabilities, we find that any equilibrium position of items subject to such forces is also 
unstable if the permittivities of all objects are higher or lower than that of the enveloping medium; 
the former being the generic case for ordinary materials in vacuum. 



Earnshaw's theorem [1] states that a collection of 
charges cannot be held in stable equilibrium solely by 
electrostatic forces. The charges can attract or repel, but 
cannot be stably levitated. While the stability of matter 
(due to quantum phenomena) , and dramatic demonstra- 
tions of levitating frogs [5], are vivid reminders of the 
caveats to this theorem, it remains a powerful indicator 
of the constraints to stability in electrostatics. There 
is much current interest in forces induced by fluctuating 
charges (e.g., mobile ions in solution), or fluctuating elec- 
tromagnetic (EM) fields (e.g., the Casimir force between 
metal plates). The former (due to thermal fluctuations) 
may lead to unexpected phenomena such as attraction 
of like-charged macroions, and is thought to be relevant 
to interactions of biological molecules. The latter (due 
mainly to quantum fluctuations) is important to the at- 
traction (and stiction) of components of microelectrome- 
chanical (MEM) devices. Here, we extend Earnshaw's 
theorem to some fluctuation-induced forces, thus placing 
strong constraints on the possibility of obtaining stable 
equilibria, and repulsion between neutral objects. 

An extension of Earnshaw's theorem [1] to polarizable 
objects by Braunbek [3J 0] establishes that dielectric and 
paramagnetic (e > 1 and /i > 1) matter cannot be sta- 
bly levitated by electrostatic forces, while diamagnetic 
(/i < 1) matter can. This is impressively demonstrated 
by superconductors and frogs that fly freely above mag- 
nets [2] ■ If the enveloping medium is not vacuum, the cri- 
teria for stability are modified by substituting the static 
electric permittivity cm arid magnetic permeability \im 
of the medium in place of the vacuum value of 1 in the 
respective inequalities. In fact, if the medium itself has a 
dielectric constant higher than the objects (e < £a/), sta- 
ble levitation is possible, as demonstrated for bubbles in 
liquids (sec Rcf . [5 , and references therein) . For dynamic 
fields the restrictions of electrostatics do not apply; for 
example, lasers can lift and hold dielectric beads with 
index of refraction n = y / e/i > 1 [6j . In addition to the 
force which keeps the bead in the center of the laser beam, 
there is radiation pressure which pushes the bead along 



the direction of the Poynting vector. Ashkin and Gordon 
have proved that no arrangement of lasers can stably lev- 
itate an object just based on radiation pressure [7]. 




FIG. 1: (Color online) Left: Each object contains mobile 
and stationary charges, which interact with charges in other 
objects according to Coulomb's law. They also interact with 
charges in the same object in an arbitrary manner and may 
be subject to an object-centered one-body potential. The 
medium has static permittivity cm- The stability of the posi- 
tion of object A is probed by displacing it infinitesimally by 
vector d. Right: The Casimir energy is considered for objects 
with electric permittivity ei(cu,x) and magnetic permeabil- 
ity fj,i(cL!,x.), embedded in a medium with uniform, isotropic, 
cm(u) and fiM(u)- To study the stability of object A, the 
rest of the objects are grouped in the combined entity R. 

We first obtain a simple extension of Earnshaw's the- 
orem to objects containing fixed and mobile charges in- 
teracting via Coulomb forces. This model, depicted in 
Fig. [I] (left), is a classical analogue of the electrody- 
namic Casimir effect, where quantum fluctuations pro- 
duce spontaneous multipole moments and currents. Does 
Earnshaw's theorem, which states that arrangements of 
fixed charges is not stable to perturbations, also apply 
to objects containing mobile charges? The free energy is 
obtained, via the partition function, by integrating the 
positions {x/} of the mobile charges {<?/} over the vol- 
umes {Vj} of the corresponding objects {J}, as 

F = -r 1 ln/ n^ e ^ ({x/}) ■ d) 
■W ev., Ij 
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where [3 — l/(fcsT). (The kinetic energy, YliPi f^ m i: is 
easily integrated out of the partition function, and since 
it is independent of the arrangement of the objects, does 
not contribute to the force.) Charges q{ and qj on differ- 
ent objects I and J interact via the Coulomb potential 
qjqjG M {x-l ,x/), where G A /(x,x') = (47re fl/ |x - x'|) _1 
is the electrostatic Green's function for a medium with 
permittivity cm, satisfying — cmV 2 Gm(x, x') = <5(x— x'). 
The Hamiltonian, neglecting the kinetic energies, is thus, 

H = E E x/) + J2 Uj ({xf }) , (2) 

KJ i,j J 



where the summation runs over i G /, j G J, and 
H , -s. 2 , . . . i represents the interactions among the 



charges and with the 'container' J. By choosing a high 
energetic cost of displacement some charges can be held 
stationary with respect to the container. 

We can determine the stability of the objects' positions 
without explicit calculation of the partition function in 
Eq. ([IJ by examining the change in free energy under 
an infinitesimal shift in the position of one object, while 
the others are held fixed. Under a translation of object 
A by d, the charges gj 4 , . . . , in A shift to positions 
x^+d, . . . , x^+d. The interaction potential Ua ({x^ 4 }) , 
however, remains unchanged since the entire container A 
is moved, that is, all vector quantities in Ua ({x^ 4 }), in- 
cluding those specifying one-body potentials, are shifted 
by +d. On the other hand, the Coulomb interaction be- 
tween a charge q^ in object A and another charge qj in 
another object J is modified to q^qj Gm (xf + d, x^) . It 
is essential that the different objects do not touch to per- 
mit the infinitesimal translation of object A. The force 
on A is given by — 'VdF. The position of object A is 
unstable if V d F < and possibly stable if V^i* 1 > 0. 

The Laplacian of the free energy is given by 



P ((V d H) 2 )-(V d H) 



(3) 



where angular brackets denote thermal averages. The 
term in square brackets equals ((Vd# — (Vd-ff)) 2 ), 
which is nonnegative and makes a destabilizing contribu- 
tion. The Laplacian in the first average only acts on the 
Green's functions in Eq. ^ which describe the interac- 
tions of the charges in object A with charges in the other 



objects, e.g., V d GA/(x^ + d, x^ )| d _ ; the (^-functions re- 
sulting from such operations, e.g., <5(x^ — x^), are always 
zero since the two charges lie in different volumes. Of 
course, the Laplacian of all other terms is zero since they 
do not depend on d. Thus, the result V d F < follows 
principally from the vanishing of V d i? for any configura- 
tion of charges within an object (as in the zero temper- 
ature Earnshaw case), and thermal fluctuations (second 
term) only enhance instability. 

Upon requiring charge neutrality for objects in the pre- 
vious model, the fluctuating mobile charges generate fluc- 



tuating electric dipoles and higher multipoles on the ob- 
jects. Even at zero temperature, quantum charge and 
current fluctuations exists, generating Casimir forces. 

Next, we proceed to this quantum mechanical case by 
considering the stability of neutral objects interacting via 
the Casimir force, which can be attributed to either fluc- 
tuations of the EM field in the medium between the ob- 
jects, or alternatively to fluctuating multipoles residing 
on them. Explicit calculations for simple geometries in- 
dicate that the direction of the force can be predicted 
based on the relative permittivities, and permeabilities, 
of the objects and the medium. Separating materials 
into two groups, with (i) permittivity higher than the 
medium or permeability lower than the medium (e > £m 
and /j, < Hm), or (ii) the other way around (e < €m 
and \i > fJ-Ai), Casimir forces are found to be attrac- 
tive between members of the same group and repulsive 
for different types. (While this has been shown in sev- 
eral examples, e.g. in Refs. U EE HU H2], a theo- 
rem regarding the sign of the force only exists for mirror 
symmetric arrangements of objects HI]-) Since or- 
dinary materials have permittivity higher than air and 
permeability very close to one, this effect causes objects 
to stick to one another. (The above statements will be 
made precise shortly.) Particularly for nanomachines this 
is detrimental as the Casimir force increases rapidly with 
decreasing separation. This has motivated research into 
reversing the force; for example, a recent experiment |15j 
shows that, in accord with the above rules, a dielectric 
medium can lead to repulsion. But the sign of the force is 
largely a matter of perspective, since attractive forces can 
be easily arranged to produce repulsion along a specific 
direction, e.g., as in Ref. [IB] . Instead, we focus on the 
question of stability which is more relevant to the design 
and development of MEMs and levitating devices. We 
find that interactions between objects within the same 
class of material cannot produce stable configurations. 
Thus, the force equilibria, as for example in Ref. [T7], are 
unstable since they are produced by geometric arrange- 
ments of infinite permittivity objects in vacuum. 

Recent theoretical advances have led to new tech- 
niques, based on scattering theory, for efficiently comput- 
ing the Casimir force (see, Ref. [TH] for a detailed deriva- 
tion, and a partial review of precursors [T31 123 121] ) ■ The 
exact Casimir energy of an arbitrary number of objects 
with linear EM response, as described in the caption of 
Fig. [l] (right), is given by (see Eq. (V.16) in Ref. [Tg]), 



£ 
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2n 



dn trlnT" 1 !. 



(4) 



where the operator [T 1 (zck, x,x')] equals 



[T_ a 1 (ick, Xi, x'j)] [<G(icK, xx, x^)] 
[G(icK, x 2 , x'x)] [T^-(icK, x 2 , x' 2 )] 



(5) 
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and Too is the inverse of T _1 with G set to zero. The 
square brackets "[ ]" denote the entire matrix or sub- 
matrix with rows indicated by x and columns by x'. 
To obtain the free energy F at finite temperature, in 
place of the ground state energy £ , f ~ is replaced 
by the sum ^ J2' K >o over Matsubara 'wavenumbers,' 
K n = 2irnkT / he, the primed summation implying that 
the «o — mode is weighted by 1/2. The operator 
)] has indices in position space. Each spa- 
tial index is limited to lie inside the objects A, B, ■ ■ ■ . 
For both indices x and x' in the same object A the 
operator is just the inverse T operator of that object, 
[T^ 1 (jck, x, x')]. For indices on different objects, x in 
A and x' in B, it equals the electromagnetic Green's 
function operator [G(ick, x, x')] for an isotropic, homo- 
geneous medium. [25] As shown in Ref. [H], after a few 
manipulations, the operators Tj and G turn into the on- 
shell scattering amplitude matrix, Fj, of object J and 
the translation matrix X, which converts wave functions 
between the origins of different objects. While practical 
computations require evaluation of the matrices in a par- 
ticular basis, the basis independent operators Tj and G 
are better suited to our general discussion here. 

To investigate the stability of object A, we group the 
'rest' of the objects into a single entity R. So, T consists 
of 2 x 2 blocks, and the integrand in Eq. Q reduces to 
tr In (I — T^GT/jG). Merging the components of R poses 
no conceptual difficulty given that the operators are ex- 
pressed in a position basis, while an actual computation 
of the force between A and R would remain a daunting 
task. If object A is moved infinitesimally by vector d, 
the Laplacian of the energy is given by 
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2T A VGT K (VG) J ftJ^ 
2T A VGT ii G I _ T J JTHG 

(TUVGT/jG + T A GT R (VG) 1 



(6) 

(7) 
(8) 



-TaGT b G 



After displacement of object A, the Green's function mul- 
tiplied by T A on the left and T R on the right (T A GT R ) 
becomes G(zck,x + d,x'), while that multiplied by 
on the left and on the right (Tr&Ta) becomes 
+ d). The two are related by transposi- 
tion, and indicated by VG(zck,x,x') = V,jG(icK,x + 
d,x')|d=o and (VG(ick, x, x')) = VdG(icK,x,x' + 
d)|d=0 in the above equation. In the first line 
we have substituted n 2 M {icK)K 2 G for V 2 G that differ 
only by derivatives of <5-functions, which vanish since 
G(jck,x,x') is evaluated with x in one object and x' 
in another. In expressions not containing inverses of T- 
operators, we can extend the domain of all operators to 
the entire space: Tj(zck,x,x') = if x or x' are not on 
object J and thus operator multiplication is unchanged. 
To determine the signs of the various terms to 



V d £| d0 , we perform an analysis similar to Ref. [IB] . 
However, we do not investigate convergence issues 
and treat the operators like matrices from the start. 
This means that the necessary criteria (smoothness, 
boundedness, compact support, etc.) are assumed 
to be fulfilled in realistic situations, as dealt with 
in Ref. [13]. The operators Tj and G are real and 
symmetric. An operator is positive (negative) semidef- 
inite if all its eigenvalues are greater than or equal 
to zero (smaller than or equal to zero). It is easy to 
verify that G is a positive semidefinite operator, since 
it is diagonal in momentum space, with G(icK,k) = 



k®k 



/(fc 2 



,{ick)k 2 



If 



f (ic K )K 2 J I V ' ,L M\ 

a real and symmetric matrix, it is positive semidef- 
inite if and only if there exists a matrix B such 
that M = B T B. Let us assume that T A and 
are each either positive or negative semidefinite, 
indicated by s A — ±1 and s R — ±1. (We shall 
shortly show how the sign of Tj can be obtained 
from the object's permittivity and permeability.) The 
eigenval ues of I — T^GT/jG, which equal those of 
I — s A \/s A T A GT R Gy/ s a T a , are strictly positive, since 
the energy is real. (The above expression appears in 
the integrand of Eq. |4]) if there are only two objects.) 
Under the trace we always encounter the combination 
(I — T^GTflG) -1 Ta, which, taking advantage of its 
symmetries and definite sign, can be written as s A M T M, 

where B = (l - s A y/s A T A GT R G \/s A T A ^ 1/2 \J s a Ta- 
The first term, line (|6]), can now be rearranged 
as trs A B T BGT fl G = s A s R tr [(BGE) (K T GB T )] 
by setting = s B RK T and its sign is s A s R . 

In the same way line 



as 



s^s^tr 



Q can be recast 
, and its sign is thus 



Lastly, the term in line ([8|) can be 

T ttdt\ 2 gince 



(BVGM) 
also set by s A s R 

rewritten as (fiVGT K GB T + BGT fl (V ; 
this is the square of a symmetric matrix, its eigenvalues 
are greater than or equal to zero, irrespective of the 
signs of T A and Tr. Overall, the Laplacian of the energy 
is smaller than or equal to zero as long as s A s R > 0. 

What determines the sign of Tjl It can be related to 
the electrodynamic 'potential' Vj, discussed in the next 
paragraph, by Tj = V,// (I + GV, 7 ) [18J. It is then pos- 
itive or negative semidefinite depending on the sign s J 



of V./, since Tj = s j t/s j Yj 



r 



l+s J y / s J Y J Gy/s 



j- 



The denominator I + is positive 

semidefinite, even if s J — —1, as its eigenvalues are the 

-l 



same as V^l.^ - ^ + Vj)vG; the term in the parantheses 
is just the (nonnegative) Hamiltonian of the field and the 
object J, G _1 +V,/ = V x ^ _1 (icK,x) V x +1 K 2 e(icn, x). 
Here, we have used Vj as given below, and e(icK, x) and 
n{icn, x) are the reponse functions defined everywhere in 
space, either of object J or of the medium, depending on 
the point x. 
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The analysis so far applies to each imaginary fre- 
quency icK. As long as the signs of and T# are 
the same over the dominant frequencies in the inte- 
gral (or the sum) in Eq. ffl, V\£ | d=0 oc —s A s R — 
(positive term). We are leftto find the sign of the po- 
tential Vj(icK,x) = I n 2 (ej(icn, x) — e^(icK)) + V x 
(icn, x) — n^(iK,)) Vx of the object A, and the 
compound object R [22] • The sign is determined by 
the relative permittivities and permeabilities of the ob- 
jects and the medium: If €j(zck,x) > ej^iicK) and 
/zj(«ck,x) < fjLMiicn) hold for all x in object J, the po- 
tential Y j is positive. If the opposite inequalities are 
true, Yj is negative. The curl operators surrounding the 
magnetic permeability do not influence the sign, as in 
computing an inner product with Yj they act symmet- 
rically on both sides. For vacuum 6m = A 1 a/ = 1, and 
material response functions e(icK, x) and /i(icK, x) are 
analytical continuations of the permittivity and perme- 
ability for real frequencies [22] - While c(zck,x) > 1 for 
positive k, there are no restrictions other than positivity 
on /z(icK,x). (For non-local and non-isotropic response, 
various inequalities must be generalized to the tensorial 
operators e (icn,x, x') and ^(ick, x, x').) 

In summary, if all objects fall into one of the two classes 
described earlier, i) cj/cm > 1 and (aj/{1m < 1 (positive 
Yj and Tj), or ii) cj/cm < 1 and Hj/hm > 1 with (neg- 
ative Yj and Tj), none of the objects levitates stably. 
(Under these conditions parallel slabs attract.) In vac- 
uum, 6M(icK) = ^m(«ck) = 1; since e(icK, x) > 1 and 
the magnetic response of ordinary materials is typically 
negligible [52~] , one concludes that stable equilibria of the 
Casimir force do not exist. If objects A and R, however, 
belong to different categories under which conditions the 
parallel plate force is repulsive-, then the terms under the 
trace in lines ^ and ^ are negative. The positive term 
in line ^ is typically smaller than the first two, as it 
involves higher powers of T and G. In this case stable 
equilibrium is possible, as demonstrated recently for a 
small inclusion within a dielectric filled cavity [23]. For 
the remaining two combinations of inequalities involving 
ej jtM and fj,j I fiM the sign of Yj cannot be determined 
a priori. But for realistic distances between objects and 
the corresponding frequency ranges, the magnetic sus- 
ceptibility is negligible for ordinary materials, and the 
inequalities involving [i can be ignored. 

For levitation in vacuum one would need a strong para- 
magnet with negligible dielectric response. Metamateri- 
als, incorporating arrays of micro-engineered circuity can 
display strong magnetic susceptibility at certain frequen- 
cies, and have been discussed as candidates for Casimir 
repulsion across vacuum j24j [25] . Given that such mate- 
rials are fabricated out of ordinary metals and dielectrics 
with well-behaved e(icK, x) and \x{ick, x) « 1 at short 
scales, they should be subject to the constraints above, 
ruling out the possibility of stable levitation. Repulsion, 
in particular, is prohibited for such materials if one of the 



objects is an infinite slab with translational symmetry; as 
the energy as a function of separation from the slab would 
then have d 2 £ > at some point since the force has to 
vanish at infinite separation. The above conclusions also 
apply to a massless scalar field, for which any arrange- 
ment of Dirichlet boundaries (corresponding to infinite 
permittivity here) is unstable to perturbations. 
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